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ABSTRACT 

Aims. We aim to detect methylamine, CH 3 NLL, in a variety of hot cores and use it as a test for the importance of photon-induced 
chemistry in ice mantles and mobility of radicals. Specifically, CH 3 NH 2 cannot be formed from atom addition to CO whereas other 
NH 2 -containing molecules such as formamide, NLLCHO, can. 

Methods. Submillimeter spectra of several massive hot core regions were taken with the James Clerk Maxwell Telescope. Abundances 
are determined with the rotational diagram method where possible. 

Results. Methylamine is not detected, giving upper limit column densities between 1.9 - 6.4 x 10 16 cnL 2 for source sizes corre¬ 
sponding to the 100 K envelope radius. Combined with previously obtained JCMT data analyzed in the same way, abundance ratios of 
CH 3 NH 2 , NLLCHO and CH 3 CN with respect to each other and to CH 3 OH are determined. These ratios are compared with Sagittarius 
B2 observations, where all species are detected, and to hot core models. 

Conclusions. The observed ratios suggest that both methylamine and formamide are overproduced by up to an order of magnitude in 
hot core models. Acetonitrile is however underproduced. The proposed chemical schemes leading to these molecules are discussed 
and reactions that need further laboratory studies are identified. The upper limits obtained in this paper can be used to guide future 
observations, especially with ALMA. 
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1. Introduction 


Complex organic molecules are thought to be formed primar¬ 


ily on dust grains in dense cores, see reviews by Herbst & van| 
Dishoeck] ( |2009| ) and |Caselli & Ceccarelli| ( |2012| >. Before the 
onset of star formation, the atomic and molecular reservoir is 
contained in large dark clouds. Due to the high densities (> 10 4 
ctrT 3 ) and low temperatures (10 K) reached in these environ¬ 
ments, gas-phase species will freeze out on sub-micron sized 
grains forming ice mantles on timescales shorter than the life¬ 
time of the cloud. It is here that atoms and molecules can poten¬ 
tially react with each other to form the zeroth order ice species 
like ammonia, methane, water and methanol. UV radiation inter¬ 
acts with these ice mantles by dissociating molecules to produce 
radicals and by photodesorbing species back to the gas phase. If 
these radicals are sufficiently mobile, they can find each other on 
the grain and react to form even more complex first generation 
(organic) species ( Garrod & Herbst|2006 1 . However, it is not en¬ 
tirely clear if UV radiation is essential to form these complex 
molecules or whether they can also be formed just by thermal 
processing and atom bombardment of solid CO with C, N and O 
atoms fTielens & Charnley|1997| >. 

In this context methylamine, CH3NH2, is a particularly 
interesting molecule, since its formation is hypothesised by 
|Garrod et al.| ( 2008) to be completely dependent on radicals 
produced by UV photons, and is one of the few molecules that 
can definitely not be produced in the routes starting from solid 
CO: 


ch 4 + hv -^ch;;+h 

NH 3 + hv -^NH^+H 

ch;+nh; - 4 CH 3 NH 2 

These radicals can form in the ice mantles in the dark 
cloud or in the protostellar phase through cosmic-ray induced 
photons and/or UV photons from the protostar. After gravita¬ 
tional collapse of the cloud and formation of a protostar, the 
dust around it will start to warm up. The increased temperature 
will cause the radicals to become mobile on the grains and react 
with each other, forming methylamine. Further heating will 
evaporate the formed methylamine from the grain and raise its 
gas-phase abundance. 

Another interesting amine-containing molecule is for¬ 
mamide, NH 2 CHO. This is so far the most abundantly observed 
amine-containing molecule (Halfen et al. 2011; Bisschop et al. 
2007 e.g.,), making it an interesting molecule to compare with 
other amines like methylamine. In contrast with CH 3 NH 2 , this 
molecule can possibly be produced by reactions of H and N 
with solid CO. The comparison of the abundances of these two 
species could potentially give more information about the rela¬ 
tive importance of UV-induced versus thermal grain surface re¬ 
actions. 

Hot cores are particularly well-suited to study methylamine. 
These high mass star-forming regions reach high temperatures 
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between 100 to 300 K and are known for their rich com- 


plex organic chemistry (Walmsley 

1992[ van Dishoeck & Blake 

1998 

Tielens & Charnley 1997; E 

irenfreund & Charnley 2000 

Case 

li & Ceccarelli 2012). The ice covered grains move inwards 


to the protostar and will heat up. When sufficient temperatures 
are reached, molecules will start to desorb depending on their re¬ 
spective binding energies. Less abundant molecules mixed with 
water ice will desorb together with water around 100 K. 

Previous detections of methylamine have all been made 


toward the galactic center. |Kaifu et al.| (1974|> first detected 


CH 3 NH 2 in Sagittarius B2 and Orion A. Later that same year 
Fourikis et al. (1974|) reported the detection of methylamine in 


the same sources, but with a different telescope. Much more sen¬ 
sitive surveys by Turner (|1991|l,|Nummelin et al. ( 2000|i,|Halfen 
|et al.| ( |20131 >, |Be loche et al.| d201 3[ > and [Neill et al.| ( |2014[ ) also 


all detected methylamine lines toward SgrB2, with typical in¬ 
ferred abundance ratios with respect to NH 2 CHO between 0.5 to 
3. No detections of methylamine have been reported in sensitive 
surveys with modem detectors toward Orion, however (Blake 


et al. 1 1987 1 |Turner| 199 1| [Sutton et al.|1995]|Schilke et al.|l997[ 


Crockett et al.|2014|>. 

To study the importance of UV processing of ice-covered 
dust grains, we present the results of searches for methylamine in 
a number of hot cores (see Table [1). These r esults are combined 
and compared with data from Bisschop et al. (2007) and Isokoski 
|et al.| ( |201 3| >, which were taken toward the same hot cores with 
the same telescope and analysis method and include detections 
of NH 2 CHO and other nitrogen-containing species. In Section 
[ 2 ] the observational details are given, followed by the analysis 
method in Section^ Section [4] summarizes all the results of our 
analysis and these are discussed in Section[5] Finally conclusions 
are drawn in Section[ 6 ] 


2. Observations 

Observations were performed with the James Clerk Maxwell 
Telescope (JCMT) Q on the sources listed in Table [T] between 
July 2010 and August 2011. The sources were selected based 
on their particularly rich chemistry, being isolated, having nar¬ 
row line widths to prevent line confusion and on their relatively 
nearby distance ( Bisschop et al.|2007)|Fontani et al.|2007||Rath-| 


|bome et al.|2008;; Isokoski et al.[2013|) . 


Nummelin et al. (1998) detected methylamine emission lines 


between 218 to 263 GHz toward Sgr B2N. Therefore the RxA3 
front-end double side band receiver, functioning between 210 to 
276 GHz, was chosen to observe the hot cores. The 250 and 1000 
MHz wide back-end ACSIS configurations were used. A number 
of methylamine transitions covering a range of excitation ener¬ 
gies were selected in this frequency range based on high Einstein 
A coefficients and lack of line confusion (Table[2|. However, not 
all transitions were observed for all sources. The 235735 MHz 
transition was only recorded for W 3 (H 20 ) and the 260293 MHz 
transition only toward W 3 (H 20 ) and NGC 7538 IRS1. 

Because double side band spectra were obtained, our spectra 
contain transitions from two different frequency regimes super¬ 
posed. To disentangle lines from the two side bands, each source 
was observed twice with an 8 MHz shift in the local oscillator 
setting between the two observations. This allows each transi¬ 
tion to be uniquely assigned to either of the two side bands. 


a The James Clerk Maxwell Telescope is operated by the Joint Astron¬ 
omy Centre on behalf of the Science and Technology Facilities Council 
of the United Kingdom, the National Research Council of Canada, and 
(until 31 March 2013) the Netherlands Organisation for Scientific Re¬ 
search. 
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Table 2. Methylamine transitions observed in this study" 


Transition 

Freq 

(MHz) 

Eup 

(K) 

A 

(s~‘) 

&up 

4 2 —> 4/ 

229310.298 

36.9 

1.32E-05 

108 

7 2 —> 7\ 

229452.603 

75.5 

5.88E-06 

60 

82 -> 8 i c 

235735.037 

92.8 

6.13E-05 

204 

62 —* 61 6 

236408.788 

60.8 

5.94E-05 

52 

2 2 —> 2i b 

237143.530 

22.0 

3.82E-05 

60 

IO 2 -» 10 / 

260293.536 

132.7 

2.26E-05 

52 


a Data from JPL database for molecular spectroscopy. 
b Transition observed in all sources. 

£ Only observed in W3(H 2 0). 
d Only observed in W3(H 2 0) and NGC 7538 IRS1. 


In the 230 GHz band, the JCMT has a beam size ( 0 B ) of 20- 
21". Spectra were scaled from the antenna temperature scale, 
T*, to main beam temperature, 7 mb, by using the main beam 
efficiency of 0.69 at 230 GHz. Integration times were such that 
Trms is generally better than 10 mK for data binned to 1.3 km 
s 1 velocity bins. Noise levels were improved by adding the 
shifted spectra together in a narrow frequency region around the 
CH3NH2 lines, effectively doubling the integration time. 


3. Data analysis 

To analyse the data, exactly the same method as described by 
|Bisschop et al.| ( |2007| ) and |Isokoski et al.| ( |2013[ > was used. It 
will be shortly reiterated here. The hot core spectra corrected 
for source velocity were analysed with the "Weeds" extension 
(Maret et al. |2011| > of the Continuum and Line Analysis 
Single-dish Software (CLASS0 coupled with the Jet Propul¬ 
sion Laboratory (JPlfl database for molecular spectroscopy 
( |Pickett et al.||1998] ). Focus was on identifying the transitions 
of methylamine listed in Table [2] but other lines in the spectra 
were measured as well (see Table |A.1| Appendix). After each 
positive identification the integrated main-beam temperature, 
f T MB dV, was determined by gaussian fitting of the line. From 
the integrated main-beam intensity the column density N up 
and thus the beam -averaged total column density At could 
be determined, assuming Local Thermodynamic Equilibrium 
(LTE) at a single excitation temperature T ml : 


3 kjT MB dV = Aup = N t , ||[|/7 . oi 
8n 3 vp 2 S g up Q(T mt ) 

where g up is the level degeneracy, k the Boltzmann constant, v 
the transition frequency, p the dipole moment and S the line 
strength. Q(T I0t ) is the rotational partition function and E up is 
the upper state energy in Kelvin. 

In case of a non-detection, 3<x upper limits were determined 
from the Root Mean Square (RMS) of the base line of the 
spectra in combination with the velocity resolution 6v and line 
width AV: 


b http://www.iram.fr/IRAMFR/GILDAS 
c http://spec.jpl.nasa.gov 
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Table 1. Source list and source parameters 


Source 

RA 

J2000 

Dec 

J2000 

AU 

AU 

L a 

Lq 

d a 

(kpc) 

Vlsr 
( km s^ 1 ) 

AV 
(km s^ 1 ) 

6v 

(km s" 1 ) 

RMS 

(mK) 

AFGF 2591 

20:29:24.60 

+40:11:18.9 

1800 

21000 

2.0E+04 

1.0 

-5.5 

4.0 

1.28 

10 

G24.78 

18:36:12.60 

-07:12:11.0 

13000 

162000 

7.9E+05 

1.1 

111.0 

6.3 

1.28 

9 

G31.41+0.31 

18:47:34.33 

-01:12:46.5 

7840 

166000 

2.6E+05 

7.9 

98.7 

7.3 

1.28 

7 

G75.78 

20:21:44.10 

+37:26:40.0 

5600 

86100 

1.9E+05 

4.1 

-0.04 

5.6 

1.28 

9 

IRAS 18089-1732 

18:11:51.40 

-17:31:28.5 

2750 

49000 

3.2E+04 

2.3 

33.8 

4.5 

1.28 

9 

IRAS 20216+4104 

20:14:26.40 

+41:13:32.5 

1753 

34400 

1.3E+04 

1.6 

-3.8 

6.0 

1.28 

10 

NGC 7538 IRS 1 

23:13:45.40 

+61:28:12.0 

4900 

58800 

1.3E+05 

2.8 

-57.4 

4.0 

1.28 

10 

W3HUO) 

02:27:04.60 

+61:52:26.0 

2400 

42000 

2.0E+04 

2.0 

-46.4 

5.0 

1.28 

11 

W33A 

18:14:38.90 

-17:52:04.0 

4500 

84000 

1.0E+05 

4.0 

37.5 

4.9 

1.28 

11 


a Data taken from Bisschop et al. (|2007 1 and Isokoski et al. (2013[>. 


cr = 1.2 VcivAU ■ RMS 


( 2 ) 


Ak is estimated from other transitions (see Table[I]i in the spec¬ 
tra, for example from the nearby FECS 1\ —> 6 | transition, and 
assumed to be the same for all transitions in the spectral range. 
A telescope flux calibration error of 20% is taken into account in 
the 1.2 factor. The 3cr value is then used in the same way as the 
main-beam intensity of detected lines to obtain the upper limit 
on the total column density through Eq. [Tj 

Since no rotational temperature can be determined for a non¬ 
detection, this has to be estimated. In the models of |Garrod et al.| 
( 20081 the peak abundance temperatures for methylamine range 
from 117 to 124 K depending on the model used. Oberg (2009 


determined that methylamine forms in CH 4 /NH 3 UV irradiation 
experiments and sublimates at 120 K. There is a small differ¬ 
ence between laboratory and hot core desorption temperatures, 
because of the pressure difference between the two. Also, if 
CH 3 NH 2 is embedded in water ice the desorption temperature 
will probably be limited to roughly 100 K, when water desorbs 
in space. Therefore T rot is assumed to be 120 K when methy¬ 
lamine lines could not be identified, but the effects of lower and 
higher rotation temperatures are explored as well. 

Correction for beam dilution is done in the same way as |Biss-| 
chop et al. ( 2007] ): 


4bf - 


fig 


e 2 s +e 2 B 


0 ) 


resulting in the source -averaged column density: 


AT ^T 

Ns - - 

r/BF 


( 4 ) 


The beam diameter 6b is set at 21". For the source diame¬ 
ter, 6s , values have been taken from Bisschop et al. ( |2007| > and 
Isokoski et al.|( 2013 ! and constitute the area where the tempera¬ 


ture is 100 K or higher and hot gas-phase molecules are present. 
Both beam and source diameters are listed in AU in TableQ] Us¬ 
ing the CASSIS line analysis software]^ it was verified that the 
source-averaged column densities are still small enough that the 
observed lines are optically thin. 


d CASSIS has been developed by IRAP-UPS/CNRS 
(http://cassis.irap.omp.eu I. 


4. Results and comparison with astrochemical 
models 


4.1. CH 3 NH 2 limits 


Figure [Tjpresents examples of spectra obtained for our sources, 
whereas Figure |A.2| in the Appendix shows the 2 2 - 2i line 
in all sources. In general, no transitions of CH 3 NH 2 are de¬ 
tected. Only one possible methylamine transition is identified in 
G31.41+0.31 coincident with the 62 —> 61 line at 236408 MHz, 
with an integrated intensity of 0.44 K kins' 1 . Following the pro¬ 
cedure summarized in Section]^] a column density of 3.4 x 10 17 
cm ~ 2 is inferred from this line assuming T [ot = 120 K. How¬ 
ever, modelling of the spectrum shows that the other targeted 
CH 3 NH 2 lines, 42 —> 4i and 22 —> 2i, should have comparable 
or even higher intensities if this identification is correct (Figure 
A. 1 1 . The 82 —» 81 line should be readily detected but was not 


observed toward G31.41+0.31. This makes it unlikely that the 
detected feature belongs to methylamine, since we would expect 
to see at least two other CH3NH2 transitions in our spectrum. 


In Figure [2] upper limit column densities of the six inves¬ 
tigated transitions of methylamine are plotted versus rotational 
temperature taking a typical 3<x = 0.100 K kms *. At 120 K the 
82 —> 81 , 235735 MHz transition gives the lowest limits on the 
column densities, see Figure [2] However, since this particular 
transition was only included in the observations for one source, 
the second most sensitive transition at 120 K, 2 2 —> 2i, will be 
used (see Figure [A~2] in the Appendix for a blow-up of this par¬ 
ticular spectral region in all investigated sources). All follow¬ 
ing molecular ratios are based on CH3NH2 column densities ob¬ 
tained from this line, assuming T rot = 120 K. The corresponding 
upper limits are presented in Table [3] 


4.2. Abundance ratio comparison 


Combined with NH 2 CHO and CH 3 OH column densities from 
Bisschop et al. (2007! and Isokoski et al.| < 2013! derived in the 
same way, abundance ratios for methylamine and formamide 
with respect to each other and to methanol are calculated. These 
ratios are listed in Table 0 Methanol is chosen as a refer¬ 
ence since it is the most readily observed complex organic 
molecule. Its disadvantage is that some of the transitions have 
high optical depth and that a cold component may be present 
(Isokoski et al.|2013|>, but this is circumvented by only taking 
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Image frequency (MHz) Image frequency (MHz) 

236800 236700 236600 236500 236400 236800 236700 236600 236500 236400 



Fig. 1 . JCMT spectra of the massive hot cores G31.41+0.31, G24.75, NGC 7538 IRS1 and W3(H 2 0). The 229310 and 229452 MHz transitions 
in the lower sideband are indicated in red and that at 236408 MHz in the upper sideband in blue. In green is the baseline, obtained by fitting line 
free portions of the spectrum. In all spectra the H 2 CS 7i —> 6i transition at 236726 MHz is fitted to determine the typical linewidth in the sources, 
as listed in Table |T| 


the warm methanol column density derived from optically thin 
lines. Abundances relative to methanol rather than H 2 are pre¬ 
ferred since the H 2 column depends on extrapolation of dust 


models to smaller scales than actually observed (Bisschop et al. 


2007). Another point that needs to be taken into account is that 
the models of Garrod et al.| ( f2008[ i do show a slight overproduc¬ 
tion of CH 3 OH, which could influence the comparison between 
the ratios. Overall, the abundance ratios are estimated to be ac¬ 
curate to a factor of a few. 


It should be noted that methylamine and formamide have 
significantly different dipole moments (1.31 and 3.73 Debye re¬ 
spectively) and could therefore be excited in different ways. For¬ 
mamide has a larger critical density than methylamine, so the 
situation could arise where the critical density is not reached for 
formamide or even both molecules. The corresponding excita¬ 
tion temperatures will then be lower. In particular, the situation 
in which the critical density is not reached for formamide but is 
for methylamine, could affect the inferred ratios. As can be seen 
from Figure [2] if T mt drops from 120 to 50 K, the column den¬ 
sity drops by a factor of a few, depending on transition. If T mt 
were 120 K for methylamine but 50 K for formamide, the ob- 
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served column density of formamide would be lower than that 
listed here and thus result in a higher C H 3 N H 2 /N H 2 C HO ratio. 
We note, however, that there is no observational evidence that 
T mt is systematically lower than 100 K for formamide (Bisschop 
|et al.|2007) . 

Table [3]inc hides the observational results toward Sgr B2, the 
only source where methylamine is firmly detected, fmm |Turner 
( fTWTT i, |Belloche et ah] ( |20T3] > and |Neill et al.| d20l4] >. These 
results, obtained over the course of more than two decades, 
agree well with each other within the estimated uncertainties 
due to slightly different adopted source sizes. Nummelin et al. 
( j2Q00) > also detect methylamine in their Sgr B2 survey but find 
a surprisingly small beam filling factor and consequently very 
large column density compared with most other complex organic 
molecules. If their beam filling factor for CH 3 NH 2 is taken to 
be the same as for NH 2 CHO, the |Nummelin et ah] ( 2000) ra- 
tios are more in line with those derived by |Turner (1991 i7 [Bel- 
|loche et aT7| ( |20131 > and |Neill et al.| ( |20T4l >. The non-detections of 
methylamine toward the chemically rich and well studied Orion 
hot core imply abundance limits that are at least a factor of 5 
lower than for SgrB2 ( |Neill et al.|20 1 4| [Crockett et al.|2014| . 
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Fig. 2. Column densities for the six methylamine transitions plotted ver¬ 
sus temperature. This plot is made for a 3<r limit of 0.1 K kms -1 , as 
found for W3(H 2 0). This figure demonstrates that the 82-81 transition 
(green) gives the most sensitive limits on column density for the rele¬ 
vant range of excitation temperatures in hot cores, when observed. The 
other five transitions (4 2 —> 4 t , black; 7 2 —> 7j, gold; 10 2 —> lOj, cyan; 
2 2 —> 2 i, blue and 6 2 —> 61 , red) clearly imply higher column densi¬ 
ties. Only below 40 K does the 2 2 —> 2 t line give lower column density 
limits. 


Table [3] also contains the model results from IGarrod et al.l 
( |2008j >, who present three hot core models which differ from 
each other by their warm-up timescale from 10 to 200 K. The 
timescales for F(ast), M(edium) and S(low) are 5 x 10 4 , 2 x 10 5 
and 1 x 10 6 years, respectively, and start after the cold collapse 
phase. In the slow models more time is spent in the warm-up 
phase where radicals are mobile. Values used in this compari¬ 
son are taken from the so-called reduced ice composition, where 
cold phase methane and methanol abundances were modified 
to match observations of these ices toward W33A, NGC 7538 
IRS9 and Sgr A*, see |Gibb et aT| ( |2000| ). Another comparison 
can be made with the gas-phase abundances in protoplanetary 
disk models of Walsh et ak| ( 2014| ) which have similar or higher 
densities and temperatures as in protostellar cores. Their ratios 
range from 7.2 x 10“' to 6.5 x 1(T 2 for CH 3 NH 2 /CH 3 OH, 4.2 x 
10 ' 1 to 1.5 for CH 3 NH 2 /NH 2 CHO and 1.7 to 8.8 x 1(T 2 for the 
NH 2 CHO/CH 3 OH. These ratios are close to the predicted values 
of |Garrod et al.|p008| listed in Table[3] which may be partly due 
to using the same surface-chemistry network. 

From Table [3] several trends become apparent for our results. 
The CH 3 NH 2 /NH 2 CHO limits lie about an order of magnitude 
above model values whereas the CH 3 NH 2 /CH 3 OH limit approx¬ 
imately matches with theoretical predictions. Because the ob¬ 
served values are actually 3 <x upper limits, this suggests that 
models overproduce CH 3 NH 2 . For the sources with the most 
stringent limits, such as G31.41+0.31 and the 8 2 - 8 ; line in 
W3(H 2 0), the CH 3 NH 2 /CH 3 OH limits are comparable or even 
lower than the abundance ratios for Sgr B2. The third ratio, 
NH 2 CHO/CH 3 OH, is also found to be lower than the models 
by up to one order of magnitude. 


Close inspection of the Bisschop et al.] (2007| data shows 
that the NH 2 CHO column densities may have larger uncertain¬ 
ties than quoted in their figures and tables. We have therefore 
re-analysed all NH 2 CHO data from that paper taking larger un¬ 
certainties into account. In general, this leads to lower NH 2 CHO 
column densities. Even using the upper limits from this re¬ 
analysis as well as those from Isokoski et al.|(|2013[) (which were 


obtained with generous error bars), the NH 2 CHO„ p/ , er /CH 3 OH 
ratios are significantly lower than the models. This suggests that 
both the methylamine and formamide abundances are too high 
in the models. 

The Sgr B2 detections tend to have lower 
CH 3 NH 2 /NH 2 CHO and CH 3 NH 2 /CH 3 OH ratios than our 
upper limits and are also somewhat below the models, but 
generally do not differ more than a factor of a few. The SgrB2 
NH 2 CHO/CH 3 OH ratios are also closer to the models results, 
at least for the faster models. However, the Orion Compact 
Ridge NH 2 CHO/CH 3 OH value from |Crockett et aT[ pOT4] ) is 
comparable to that found for our sources and clearly lower than 
the models. Further observations are needed to determine to 
what extent Sgr B2 is a special case. 

To further elucidate the differences between theory and our 
and the Sgr B2 observations, an additional analysis was made 
of the CH 3 NH 2 /CH 3 CN ratio. These results are listed in Table 
[4] Acetonitrile is produced in the gas-phase, but more abun¬ 
dantly on grains: an important route to its formation is via CH* + 
CN* —> CH 3 CN, according to Garrod et al. ( 2008[ ). This would 
mean that both molecules compete for the methyl radical on the 
surface, thus relating the two molecules. 

Our observed ratios involving CH 3 CN are clearly at odds 
with the theoretical predictions. The observed CH 3 NH 2 /CH 3 CN 
ratios are in most cases an order of magnitude lower than the¬ 
ory and approach the observed ratios for Sgr B2. However, the 
observed CH 3 CN/CH 3 OH ratios are 1-2 orders of magnitude 
higher than theoretical predictions. Both these cases point to 
CH 3 CN being underproduced in the models. 

Finally, abundance ratios, with some notable exceptions, do 
not vary more than an order of magnitude between different 
sources, as also found by Bisschop et al. (2007) for other species. 


5. Discussion 

Despite a significant number of succesfully identified molecules 
(see Table [ATT] in the Appendix for examples of W3(H 2 0) and 
G31.41+0.31), only upper limits were found for methylamine in 
the various hot cores, limiting the conclusions that can be drawn. 
Nevertheless, trends are seen in our abundance ratios. The results 
suggest that theoretically predicted abundances for both methy¬ 
lamine and formamide are too high. In contrast, acetonitrile is 
found to be underproduced in the models. In the following, each 
of these species is discussed individually. 


5.1. CH 3 NH 2 


Garrod et al. ([2008 1 suggest that methylamine is primarily 


formed by grain-surface chemistry using UV to create the CH 3 
and NH 2 radicals from photodissociation of primarily CH 4 and 
NH 3 . Perhaps the amount of UV processing is overestimated 
in these models. An alternative route is hydrogen atom addi¬ 
tion to solid HCN, proposed by |Theule et al.| ( |2011) > and found 
to lead to both CH 2 NH (methanimine) and CH 3 NH 2 . |Walsh| 
|et al.| ( [20T4| find in their models that methylamine is indeed ef¬ 
ficiently formed on grains at 10 K by atom addition reactions 


to solid CH 2 NH. Burgdorf et al. (2010! have detected HCN ice 
on Triton, but so far no detection of solid HCN has been made 
in the ISM. Methanimine is actually readily observed in the 
gas-phase (|Turner|[i991 Nummelin et al.p998| 


Belloche et al. 


|2013| > so the presence of both species makes the H-atom addition 
scheme probable. However, Halfen et al. ( [2013) detect CH 2 NH 
in Sgr B2(N) at a rotational temperature of 44 K, which is dis¬ 
tinctly colder than the 159 K observed for CH 3 NH 2 , suggesting 
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Table 3. Upper limit column densities and abundance ratios for methylamine. 


Source 

Ns.CHjNHj 

cm~ 2 

CH 3 NH 2 /NH 2 CHO 

CH 3 NH 2 /CH 3 OH 

NH 2 CHO/CH 3 OH 

NH 2 CHO u; , ;) „,/CH 3 OH 

Model F 


1.1 

3.4E-02 

3.1E-02 

3.1E-02 

Model M 


1.7 

1.0E-01 

7.3E-02 

7.3E-02 

Model S 


1.3 

1.3E-01 

1.0E-01 

1.0E-01 

AFGL 2591 

<1.9E+16 

- 

- 

<3.9E-01 

- 

G24.78 

<2.4E+16 

<3.3E+01 

<8.5E-02 

2.6E-03 

9.0E-04 

G31.41+0.31 

<5.8E+16 

<2.8E+01 

<4.9E-02 

1.8E-03 

3.8E-03 

G75.78 

<3.5E+16 

<1.7E+02 

<3.1E-01 

1.8E-03 

2.6E-02 

IRAS 18089-1732 

<4.2E+16 

<5.0E+01 

<1.9E-01 

3.8E-03 

7.9E-03 

IRAS 20216+4104 

<6.4E+16 

- 

<2.2 

- 

- 

NGC 7538 IRS 1 

<2.0E+16 

<3.5E+01 

<1.8E-01 

4.8E-03 

2.1E-04 

W3(HoO) 

<5.0E+16 

<3.9E+01 

<5.0E-02 

1.3E-03 

6.4E-04 

W3(H 2 0)* 

<1.7E+16 

<1.3E+01 

<1.7E-02 

1.3E-03 

6.4E-04 

W33A 

<5.7E+16 

<2.7E+01 

<2.9E-01 

1 .1E-02 

4.6E-03 

Sgr B2 a 

1.2E+14 

5.7E-01 

1.7E-02 

1.3E-02 


Sgr B2(M) b 

4.5E+16 

3.2 

1.7E-02 

5.2E-03 


Sgr B2(N)* 

6.0E+17 

4.3E-01 

3.3E-02 

7.8E-02 


Sgr B2(N) C 

5.0E+17 

2.1 

1.0E-01 

4.8E-02 


Orion Compact Ridge' 1 ' 

- 

- 

- 

1.6E-03 



Notes. Column densities for the assumed source size and upper limit abundance ratios for methylamine, derived from the 2 2 - 2i line assuming T rot 
= 120 K. The values for NH 2 CHO and CH 3 OH where taken fr om|Bisschop et aT7| ( [2007| > and |Isokoski et al.|(|2013|>. The u pper limits of NH 2 CHO 
were determined by our own re-analysis of the Bisschop et al. ( 2 007) data and taken from the appendix of Isokoski et al. (2 013^ . 

* Column de nsity calculated for the 8 2 - 8i line. 


References. “ 


Turner 


11991 1 , beam sizes between 65" and 107", assuming no beam dilution; b Belloche et al.| l 2013 i, beam sizes between 25" and 
30", assu ming a 3" source size for (N) and 5" source size for (M); c |Neill et ah|(2014) , beam srzes between 10" and 40", assuming source size of 
2.5"; and 0 Crockett et al. (20141 beam sizes between 44" and 11" and assuming a 10" size of the Compact Ridge. 


that the two molecules may not co-exist. An alternative route 
would therefore be to form these molecules by two different gas- 
phase reaction pathways (CH*(g) + NH 3 (g) —> CH 2 NH + H and 
CH*(g) + NH 3 (g) -> CH 3 NH 2 + H), with CH being present pri¬ 
marily in the colder outer envelope and CH 3 in the warmer cen¬ 
ter. Further modeling is needed to determine whether these gas- 
phase reactions can reproduce the observed abundances quanti¬ 
tatively. 


5.2. NHsCHO 


Formamide also appears to be overproduced in the hot core 
model. Since !Garrod et al. ( |2008 1 use both gas-phase, radical 
and atom addition reactions to form formamide, it is difficult to 
pin down where the discrepancies could come from. It is known 
that NH 2 CHO is formed in CO:NH 3 mixtures after UV and elec¬ 
tron irradiation ( |Grim et al.|1989|[Demyk et al.|1998||Jones et ah] 
|2011[ l and it has also been proposed that it can form from H- and 
N-atom addition to solid CO ( |Tielens & Charnley|jl997 1 . Gas- 
phase formation from CO and NH 3 is viable as well ( |Hub bard 
|et al.||T975) , although these experiments were conducted under 
high-pressure conditions, not the low pressures applicable in the 
ISM. Further quantification of both gas-phase and solid phase 
routes through laboratory experiments is needed. Recent labora¬ 
tory experiments by Fedoseev et al. ( 2014] l do not find NH 2 CHO 
production in H- and N-atom bombardment studies of solid CO, 
consistent with a large barrier for H- addition to HNCO found 
in ab initio calculations ( Nguyen et al.|1996] >, so perhaps the ef¬ 
ficiency of this route has been overestimated in the models. An 
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alternative solution would be that the high-mass sources studied 
here have not gone through a long (pre-stellar) phase in which 
the dust temperature was low enough for CO to be frozen out 
and turned into other molecules. 


5.3. CH.CN 


The clear mismatches between theory and observations for the 
ratios involving CH 3 CN point toward an underproduction of ace¬ 
tonitrile by more than an order of magnitude in the models. As 
with formamide, gas-phase, radical and atom addition reactions 
contribute to the formation of CH 3 CN in the models, making it 
difficult to determine the cause. The main formation route in the 
models by radical addition of solid CH* and solid CN* has never 
experimentally been investigated. It would therefore be useful to 
determine if this is a viable solid state formation route and if it 
potentially has a higher efficiency than assumed. 

Alternatively it is possible that photodestruction of solid ace¬ 
tonitrile is not as efficient as assumed in the models. Gratier et ahl 
( j2013) > find high gas-phase CH 3 CN abundances in the Horse- 
head PDR, indicative of a high photodesorption rate and slow 
destruction of CH 3 CN in the ice. Bernstein e t ak] ( j2004| > indeed 
find slower photolysis of solid CH 3 CN compared with other or¬ 
ganic molecules. If such a slower photodissociation rate would 
also hold for gas-phase CH 3 CN, it would be an attractive ex¬ 
planation why the CH 3 CN rotational temperatures are gener¬ 


ally higher than those of other complex molecules (e.g., Biss 
chop et al. (|2007]> and many other hot core studies), since the 


molecule could then approach the protostar closer before being 
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Table 4. Upper limit column densities and abundance ratios for methy¬ 
lamine. 


Source 

CH 3 NH 2 /CH 3 CN 

CH 3 CN/CH 3 OH 

Model F 

5.5E+01 

6.3E-04 

Model M 

3.8E+01 

2.6E-03 

Model S 

1.5E+01 

8.6E-03 

AFGL 2591 

<5.3 

<7.5E-02 

G24.78 

<5.1E-01 

2.1E-01 

G31.41+0.31 

<3.6 

5.9E-02* 

G75.78 

<1.9E+01 

1.6E-02 

IRAS 18089-1732 

<8.9 

1.1E-02* 

IRAS 20216+4104 

<4.3E+01 

5.2E-02 

NGC 7538 IRS 1 

<2.5 

6.8E-02 

W3(H->0) 

<7.2 

7.0E-03 

W33A 

< 2.1 

1.4E-01 

Sgr B2 fl 

1.2 

1.5E-02 

Sgr B2(N) i 

3.0E-01 

1.1E-01 

Sgr B2(M) i 

2.5E-01 

6.7E-02 

Sgr B2(N) f 

5.9E+01 

1.7E-02 

Orion Compact Ridge 4 

- 

1 . 1 E -02 


Notes. Abundance ratios for methylamine. The values for Cl LOH and 
CH 3 CN were taken from Bisschop et al.| (|2007|> and |Isokoski et al.| 

*Ratio derived from optically thi n C isotope. 


References. " Turner 
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Id 199 it. 


( |2014j > and rf |Crockett et al.H2014| T 


Belloche et al. (2013 s c Neill et al. 


destroyed. However, the gas phase photoabsorption cross sec¬ 
tions of CH 3 CN are well determined and if the bulk of these ab¬ 
sorptions lead to dissociation this would result in a photodisso¬ 
ciation rate of gaseous CH 3 CN at least as fast as that of CH 3 OH 
fvan Dishoeck et al.|2006j >. 

Another important parameter for all molecules studied here 
is the mobility of radicals and neutral molecules on the sur¬ 
face assumed in the gas-grain models. For many species no 
experimental data are available on diffusion barriers, only 
theoretically-inspired guesses. Observational evidence suggests 
that at least parts of the ice mantles are segregated in CO-rich and 
CO-poor layers (Helens et al. 1 1 99~f] |Pontoppidan et al.||2008] l. 
Therefore, more knowledge of the structure of ice mantles and 
the mobility of radicals and neutral molecules as a function of 
surface temperature and in various chemical environments is 
necessary to determine if addition reactions are likely to happen 
and at which rates. 


5.4. Prospects for ALMA 

In the near future, much deeper searches for CH 3 NH 2 can be 
carried out by ALMA (Appendix |Bj). Figure [B~T| shows that the 
strongest transitions within Band 6 are mainly located between 
240 and 275 GHz and in Band 7 around 310, 340 and 355 GHz. 
In Table B. 1 1 the strongest transitions in ALMA’s Band 6 and 7 
are listed. It becomes apparent that lines covered by Band 7 are 
more intense, but at the cost of a lower line density. 

Estimates done for the W 33A source with the CASSIS line 
analysis software and the ALMA Sensitivity Calculator show 
that ALMA should be able to reach the 3cr detection limits for 
the CH 3 NH 2 lines around 236 GHz in less than 1 hour of inte¬ 


gration time, assuming the column density for methylamine of 


1 . 2 xl 0 14 cm 2 as found by 


Turner 


(19911 in a large beam and 


two orders of magnitude lower than those inferred here for a 
small source size. This estimate assumes a spectral resolution of 
0.64 kms 1 as used in our JCMT data, the number of ALMA an¬ 
tennas set to 34 (as in Cycle 2) and a synthesized beam of 1.1", 
appropriate for the W33A hot core (100 K radius). 


6. Conclusions 

We have analysed nine hot core regions in search of methy¬ 
lamine. The molecule has not been convincingly detected, so 
upper limit abundances are determined for all the sources. 
From these limits, ratios of methylamine to other molecules 
(NH 2 CHO, CH 3 OH, CH 3 CN) have been determined and com¬ 
pared with theory and Sagittarius B2 surveys. Our conclusions 
are as follows: 

1. Trends in our results indicate that both methylamine and for- 
mamide are overproduced in the models of |Garrod et al.| 
( 2008[ >. Acetonitrile is underproduced with respect to these 
models. This is especially true for the slow models. 

2. Abundance ratios do not differ more than an order of mag¬ 
nitude between various sources suggesting that the (nitro¬ 
gen) chemistry is very similar between hot cores, as has been 
found previously for other species. 

3. More (laboratory) studies are needed to clarify the formation 
pathway of methylamine and to determine differences and 
similarities with formamide, methanimine and, to a lesser 
extent, acetonitrile formation. 

4. The upper limits determined for CH 3 NH 2 here can guide 
future more sensitive observations, especially with ALMA. 
Based on the ratios found in the Sgr B2 observations it is very 
likely that ALMA will reach the detection limit for methy¬ 
lamine in the sources studied here. Particularly strong transi¬ 
tions and spectral regions to target with ALMA are given. 
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Appendix A: Methylamine and other transitions 



|T MB dV (K kms- 1 ) 
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-!-TT--1--1- * -I- -I- r--1- 

42-4i 72-7i 82-81 62-61 22 -2 i IO2-IO1 


Transition 


Fig. A.l. Bargraph plot of the integrated main-beam intensities for the six investigated methylamine transitions, f T MB dV values were calculated 
for a total column density of 3.4 x 10 17 cnT 2 as inferred toward G31.41+0.31 from the 62 —> 61 transition, assuming a rotational temperature of 
120 K. The horizontal line shows the 0.077 K kms -1 3<r value for G31.41+0.31. It is unlikely that the detection of the 236408 feature (green) in 
this hot core is methylamine, since the main-beam intensities of other CH 3 NH 2 transitions are above the 3 cr value. Particularly the 42 —> 4i and 2i 
—> 2! transition (blue) should be visible in our spectra. The remaining transitions (red) are either below detection limit or not observed toward this 


source. 
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Table A.l. All identified transitions for the sources G31.41+0.31 and W3(H 2 0), with integrated peak area listed. 


Species 

Freq 

(MHz) 

£ up 

(K) 

A 

(s- 1 ) 

Transition G31.41+0.31 W3(H 2 0) 

(Kkms- 1 ) (K kms- 1 ) 

CH 3 OCHO 

228628.876 

118.8 

1.66E-04 

I 85132 —> 175122 

- 

0.39 

CH 3 OCHO 

228651.404 

118.8 

1.66E-04 

I 85130 —> 175120 

- 

0.53 

CH 3 COCH 3 

228668.358 

85.4 

1.69E-04 

14961 —> 13852 

- 

0.24 

CH 3 OCHO 

229388.947 

217.0 

1.62E-05 

239150 —> 238160 

0.79 

- 

CH 3 OCHO 

229405.021 

110.7 

1.75E-04 

183152 —> 173142 

2.18 

0.78 

CH 3 OCHO 

229420.342 

110.7 

1.75E-04 

183150 ~> 233140 

1.86 

0.71 

CH 3 OCHO 

229504.724 

134.3 

1.18E-05 

2 O 3170 —> 194160 

1.52 

- 

CH 3 OH 

229589.056 

374.4 

2.08E-05 

1540 —> I 630 

2.58 

1.06 

CH 3 OH 

229758.756 

89.1 

4.19E-05 

810 —>7oo 

7.17 

- 

CH 3 OCHO 

236355.948 

128.0 

1.93E-04 

20318 1 — * 19317 1 

2.38 

- 

CH 3 OCHO 

236365.574 

128.0 

1.93E-04 

203180 —» 193170 

1.93 

- 

CH 3 NH 2 

236408.788 

60.8 

5.94E-05 

621 —» 610 

0.44 

- 

HCCCN 

236512.777 

153.2 

1.05E-03 

26 -> 25 

4.67 

2.27 

h 2 cs 

236726.770 

58.6 

1.91E-04 

7l 7 —» 61 6 

7.54 

5.75 

CH 3 OCHO 

236743.697 

129.6 

1.86E-04 

19515 1 —> I 8514 1 

2.21 

?« 

CH 3 OCHO 

236759.687 

129.6 

1.86E-04 

195150 ~** 185140 

1.40 

0.53 

CH 3 OCHO 

236800.589 

136.7 

1.80E-04 

19614 1 —> 18613 1 

1.21 

? 

CH 3 OCHO 

236810.314 

136.7 

1.81E-04 

196140 ~> 186130 

2.50 

1.20 

CH 3 OH 

236936.089 

260.2 

2.79E-05 

14-10 —> 1320 

2.74 

1.91 

CH 3 OCHO 

236975.844 

320.3 

2.01E-04 

22 122 4 —> 21 1214 

0.86 

0.33 

CH 3 OCHO 

236976.390 

320.3 

2.01E-04 

22022 5 —> 2 I 0210 

0.86 

0.33 

CH 3 OCHO 

236975.844 

320.3 

2.01E-04 

22 122 4 —> 21 121 4 

1.19 

0.52 

CH 3 OCHO 

236976.390 

320.3 

2.01E-04 

22 o 22 5 ~> 21o215 

1.19 

0.52 

CH 3 OCH 3 

237046.092 

31.3 

2.33E-05 

7253 ~> 61 6 3 

2.90 

0.80 

CH 3 OCH 3 

237046.106 

31.3 

2.33E-05 

7252 —>6162 

2.90 

0.80 

CH 3 OCH 3 

237048.797 

31.3 

2.32E-05 

725 1 —^ 6 l 6 1 

2.90 

0.80 

CH 3 OCH 3 

237051.495 

31.3 

2.33E-05 

7250 ~> 61 60 

2.90 

0.80 

S0 2 

237068.870 

94.0 

1.14E-04 

123 9 ~> 122 10 

1.07 

1.96 

OC 34 S 

237273.635 

119.6 

3.88E-05 

20-» 19 

1.24 

0.58 

CH 30 CHO 

237297.482 

128.0 

1.95E-04 

20218 2 —> 19217 2 

? 

2.45 

CH 30 CHO 

237309.540 

131.6 

1.98E-04 

2 1 2201 ~> 20219 1 

? 

2.45 

CH 30 CHO 

237315.082 

131.6 

1.98E-04 

2 1 220 0 ~> 20219 0 

? 

2.45 

CH 30 CHO 

237344.870 

131.6 

1.98E-04 

21 l 202 —> 20h92 

3.48 

1.19 

CH 30 CHO 

237350.386 

131.6 

1.98E-04 

21 l 200 ~> 20„9o 

3.48 

1.19 

ocs 

243218.040 

122.6 

4.18E-05 

20 -» 19 

- 

2.95 

CH 30 CH 3 

259982.561 

226.6 

7.27E-05 

2 O 5162 —> 2 O 4172 

- 

1.54 

CH 30 CH 3 

259982.596 

226.6 

7.27E-05 

2 O 5163 2 O 4173 

- 

1.54 

CH 30 CH 3 

259984.480 

226.6 

7.27E-05 

2 O 5161 —> 2 O 4171 

- 

1.54 

CH 30 CH 3 

259982.561 

226.6 

7.27E-05 

2 O 5160 —> 2 O 4170 

- 

1.54 

NH 2 CHO 

260189.848 

92.4 

1.25E-03 

12210 —» 1129 

- 

1.46 

h 13 co+ 

260255.339 

25.0 

1.33E-03 

3 -+ 2 

- 

11.76 

CH 30 CH 3 

260327.165 

208.3 

7.21E-05 

195152 —> 194162 

- 

1.27 

CH 30 CH 3 

260327.238 

208.3 

7.21E-05 

19515 3 ~» 19416 3 

- 

1.27 

CH 30 CH 3 

260329.312 

208.3 

7.21E-05 

19515 1 ~> 19416 1 

- 

1.27 

CH 30 CH 3 

260331.422 

208.3 

7.21E-05 

195150 ~> 194160 

- 

1.27 

SiO 

260518.020 

43.8 

7.21E-05 

6 —> 5 

- 

9.24 


Indicates an identified transition where it was not possible to determine the peak area. 
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Image Frequency (MHz) Image Frequency (MHz) 

237050 237100 237150 237200 237100 237150 237200 237250 



Fig. A.2. Blow-up of the spectral region around the 22 —> 2i transition at 237143 MHz of all analysed hot cores. Despite being a particularly strong 
transition, it was not observed in any of the spectra. 
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Appendix B: ALMA 

In the following figure the simulated spectrum of methylamine is shown for a column density of 1.0 X 10 15 cm~ 2 at an excitation 
temperature of 120 K. Simulations were done with CASSIS using the JPL spectroscopic database. The frequency ranges were taken 
to cover ALMA Bands 6 and 7. The resolution was set at 0.1 MHz for this spectrum. 



Spectral Flux Density (Jy) 


A&A-Methylamine_hotcores_Ligterink, Online Material p 14 



220 240 260 280 300 320 340 360 


Frequency (GHz) 


Fig. B.l. 211 to 373 GHz spectrum of methylamine, covering ALMA bands 6 and 7. A column density of 10 15 cm 2 , T rot and beam size of 1" are 
used. 

Table B.l. Methylamine transition target candidates for ALMA Band 6 and 7 


Transition 

Freq 

(MHz) 

£u P 

(K) 

A 

(s- 1 ) 

72 6 -* 7i 6 

227498.06 

75.4 

4.80E-05 

52 6 -> 5i 6 

232003.95 

47.7 

4.73E-05 

82 3 —^ 81 2 

235735.04 

92.8 

6.13E-05 

72 2 —> 7i 3 

239209.63 

75.8 

6.29E-05 

62 3 —> 61 2 

242262.02 

60.9 

6.39E-05 

52 2 -> 5i 3 

244886.90 

48.1 

6.42E-05 

80 5 -» 7i 5 

259042.46 

77.0 

5.73E-05 

62 2 ~> 61 3 

261252.89 

60.9 

7.14E-05 

80 2 7j 3 

261563.15 

76.8 

5.98E-05 

4l 5 3o 5 

264172.21 

25.9 

8.74E-05 

CO 

00 

T 

<N 

Cd 

OO 

268898.14 

92.8 

7.44E-05 

9o 5 —> 8 j 5 

299189.80 

96.1 

8.78E-05 

9o 3 —> 81 2 

301654091 

95.9 

9.11E-05 

5-t 3 -> 4o 2 

307791.75 

36.3 

1.44E-04 

14 2 2 —> 14_i 3 

310750.84 

240.0 

8.33E-05 

10o 5 —> 9i 5 

338628.32 

117.3 

1.25E-04 

10 o 2 —> 9_i 3 

341059.48 

117.1 

1.30E-04 

CO 

0 

T 

(N 

VO 

354843.73 

49.2 

2.16E-04 

615 -> 5o 5 

357440.12 

49.5 

2.16E-04 

32 4 —> 2 i 4 

370166.34 

28.4 

2.11E-04 



